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ABSTRACT 

The present work represents an attempt to understand the 'rules of behavior' of 
observed warps in the HI disks of spiral galaxies found by Briggs (1990). In contrast 
with most earlier theoretical work, the present study investigates different initial value 
problems of a warped disk in an oblate (or prolate) halo potential, and it represents the 
disk warp in terms of N independently tilted, self-gravitating, concentric rings. This 
representation gives new insight into the disk warping. A new constant of the motion 
of N tilted rings is identified (in addition to the energy). The phenomenon of phase- 
locking of the lines-of-nodes of nearby rings due to self-gravity is demonstrated. We 
consider the influence of dynamical friction due to ring motion through the halo matter 
as well as friction between gaseous rings with different vertical motions due to turbulent 
viscosity. We first consider the dynamics of one, two, and three tilted rings of different 
radii in a halo potential. We go on to develop dynamical equations for N— rings which 
are most simply expressed in terms of the complex tilt angles Qj = 9jexp(—iipj), where 
6j is the actual tilt angle and tfj the line-of-nodes angle for the j th ring (j = 1..N). We 
numerically solve the equations for Qj for four different types of initial conditions: (1) 
warp excitation by a passing satellite, (2) excitation by a sinking compact minor satellite, 
(3) warp evolution due to a tilted halo potential, and (4) warp evolution resulting from 
an initially tilted disk plane. 

Subject headings: accretion, accretion disks — instabilities — galaxies: kinematics and 
dynamics 
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1. Introduction 

The large-scale integral-sign warps of the outer re- 
gions of spiral galaxies has been a long-standing the- 
oretical puzzle. The strongest warps are observed in 
the HI disks when these extend well beyond the opti- 
cal disks (Briggs 1990). The warps are antisymmetric 
and consequently can be well-fitted by the kinematic 
tiltcd-ring model of Rogstad, Lockhart, and Wright 
(1974). For a set of 12 galaxies with high-quality 
data on the extended HI disks, Briggs (1990) found 
empirically a set of 'rules of behavior' of warps: 1. 
The warps develop with increasing radius r for r be- 
tween i?25 and RhoI (the Holmberg radius, = i?26.5)> 
where the subscript denotes the surface brightness in 
B in mag/arcsec 2 ; 2. The line-of-nodes of the warp is 
straight for r < RhoV, and 3. For r > RhoI, the line- 
of-nodes forms an open leading spiral. An illustrative 
case is shown in Figure 1 which gives a surface plot 
of the warp of the disk of the spiral galaxy M 83 ob- 
tained from the data of Briggs (1990). The rotation of 
M 83 is clockwise so that the warp is a leading spiral 
wave. Figure 2 also shows the data of Briggs (1990) 
for M 83, in panel (a) the radial (r) dependence of 
the tilt angle 6, in (b) the r-dependence of the angle 
of the line-of-nodes ip, and in (c) the dependence of 9 
on ip, which we refer to as a Briggs plot. A different 
behavior is observed for the galaxy NGC 3718 which 
shows an approximately straight line-of-nodes out to 
r w 1.75RhoI and a large tilt angle (~ 90°) at this 
radius (Briggs 1990). 

A partial listing of the theoretical works on warps 
of disk galaxies include Kahn and Woltjer (1959), 
Hunter and Toomre (1969; hereafter HT), Toomre 
(1983), Sparke and Casertano (1988), Binney 1992, 
Hunter (1994), Merritt and Sellwood (1994), and Nel- 
son and Trcmaine (1995). Since the pioneering anal- 
ysis by Hunter and Toomre (1969) and the earlier 
related work by Lynden-Bell (1965), the theoretical 
studies have focused mainly on determining the eigen- 
modes and eigenfrequencies of warped self-gravitating 
disks. The studies used the HT dynamical equa- 
tion for the vertical displacement of the disk surface 
h(r,(f),t), where <j> is the azimuthal angle. 

In contrast with these earlier works, the present 
study investigates different initial value problems of a 
warped disk in an oblate (or prolate) halo potential, 
and it develops a different representation for the disk 
warping. Instead of describing the disk in terms of 
h(r,4>,t), we develop a representation in terms of N 



independently tilted, self-gravitating, concentric rings 
which is suggested by the kinematic ring model of 
Rogstad et al. (1974). May and James (1984) men- 
tion a model of this kind, but they give no details. We 
show that the N— ring representation is equivalent to 
the HT equation for h{r,<j),t) for N — > oo. However, 
new insight is given by the TV— ring model: A new 
constant of the motion for a warped disk is found. 
The notion of the phase-locking of the line-of-nodes 
of nearby rings due to self-gravity is demonstrated. 
We propose that this explains Briggs's rule No. 2. 
Furthermore, treatment of initial value problems for 
the age of a galaxy can be done with sufficient accu- 
racy with relatively small values of N (~10 2 ). 

In §2 we first give basic parameters of the disk and 
the halo. In §2.1 we treat the case of one tilted ring 
and in §2.2 comment on the interpretation of obser- 
vations of a ring. Subsection 2.3 considers the case of 
two tilted rings. Subsection 2.4 treats the case of two 
tilted, counter-rotating rings. Subsection 2.5 treats 
the case of three corotating rings. Section 3 treats 
the general case of N tilted corotating rings. In this 
Section we consider different possible excitations or 
origins of warps - that due to a passing satellite (§3.1), 
a sinking satellite (§3.2), a tilted halo potential (§3.3), 
and an initially tilted outer disk plane (§3.4). Section 
4 discusses the continuum limit where N — ► oo. Sec- 
tion 5 gives conclusions of this work. 

2. Tilting of Disk Galaxies 

The equilibrium, unwarped galaxy is assumed to 
be axisymmetric and to consist of a thin disk of stars 
and gas and a slightly oblate or prolate halo of dark 
matter. We use an inertial cylindrical (r, <f>, z) and 
Cartesian (x, y, z) coordinate systems with the disk 
and halo equatorial planes in the z — plane. The 
total gravitational potential is written as 



$(r, z) = <5> ld + $ od + $ h , 



(1) 



where $^ is the potential due to the inner part of the 
disk (as discussed below), $ 0( j is that for the outer 
disk, and <&h is that for the halo. The particle orbits in 
the equilibrium disk are approximately circular with 
angular rotation rate O(r), where 



i a$ 

r dr 



= ^ d + n 2 od + n 2 h 



(2) 
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The vertical epicyclic frequency (squared), 



KzW " dz 2 



K ldz 



K hz ) 



(3) 



gives a measure of the restoring force in the z— direc- 
tion. 

The surface mass density of the inner (optical) disk 
is assumed to be £ = E exp(— r/r<j) with S and 
constants and Md = 2Trr 2 T, the total disk mass. For 
r 2 > r 2 , 

GM d 



2 _ GM d 

K idz — q 
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27 r 2 



(4) 



(Binney and Tremaine 1987; hereafter BT, p. 409). 
Note that^-fi 2 d = 9GM d r 2 /r 5 . 

The outer, neutral hydrogen disk can be described 
approximately by a Fermi function for the surface 
density 



-'H 



l + exp[fc!(r/i? ff -k 2 )} 



(5) 



where Rh is the radius inside of which half the neutral 
hydrogen mass is located, k\ w 3.70, and ~ 1.21. 
This dependence provides a good fit to the data of 
Brocils and Rhce (1966; their Figure 6a). The total 
neutral hydrogen mass is Mu 

The halo potential is taken to be 



5.35R 2 H Z H0 . 



1 



®h = Wo in [r 2 + r 2 + (z/qf 



(6a) 



where v Q = const, is the circular velocity at large 
distances, r a = const, is the core radius of the halo, 
and q = const, is the axial ratio of the equipotential 
surfaces with s = 1 — q their ellipticity. An oblate 
(prolate) halo corresponds to q < 1 (q > 1). We have 



K hz 



(Vo/qf 



(66) 



The halo is usually assumed to be oblate with e ~ 
0.1 — 0.2 (Nelson and Tremaine 1995). This gives 

K L ~ n h = (V~ 2 - ~ (°- 23 ~ °- 56 )^ft which is 
of interest in the following. For the radii of the warps 

observed in galaxies (r ~ 4r"d), n 2 dz — Qf d <C ft 2 ,- For 
q = 1, the halo mass within a radius r is Mh(r) w 
0.925 xlO n M (r/lOkpc)(u o /2OOkm/s) 2 [r 2 /(r 2 + r2)]. 

Perturbations of the galaxy are assumed to consist 
of small angle tilting (9 2 <C 1) of the outer disk with 



azimuthal mode number m = 1 (or —1). The halo 
and the inner disk are assumed to be unaffected by 
the perturbation. We can describe the outer disk by 
a number ./V of tilted plane circular rings. This de- 
scription is general for small tilting angles where the 
linearized equations are applicable. We do not need 
to assume that the disk is 'razor thin' but rather that 
the disk half-thickness Az is small, (Az) 2 <C r 2 . The 
vertical displacements of the disk ~ Or are in general 
much larger than Az. 

2.1. One Tilted Ring 

The tilting of the ring is completely described by 
the tilt angle 0(t), which is the angle between the up- 
ward normal to the ring and the z— axis, and the az- 
imuthal angle <p(t), which is the angle between the 
line-of-nodes (where the ring intersects the z — 
plane) and the a;— axis. The geometry is shown in 
Figure 3. The perturbation is assumed to be small 
in the respect that 2 <C 1. The position vector of a 
point on the ring is 



R = [ rcos0, rsin0, h(cj>, t) 



< <b < 2tt 



(7) 



(Cartesian components), where h(<f>,t) is the height 
of the ring above the z = plane. We have 



h(cj>, t) = r9(t)sm [<j> - tp(t)} , 
= r0 x (t)sa\<p — r8 y (t)cos(j) 



where 



9 x (t) ee 9(t)cos [<p(t)] 
is the angle of the tilt about the x— axis, and 

y (t) ee 0(i)sin W)] 



(8a) 
(86) 

(8c) 



is the angle about the axis. We have 6 = ^Jo 2 + 9 2 
and <p = ta,TT 1 (9y/9 x ). 

The angular momentum of the tilted ring is 



. , f 27T # dn 

M — R x 

Jo 2tt 



dt 



= Mnr 2 z + 5L 



(9) 



where = xSL x + ySL y , M is the ring mass, and f2 
is its angular rotation rate which is unaffected by the 
perturbation to first order in 9. We have 



dR 

dt 



ds 



ds &R 
dt + ~dt 
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tftr 



<9R 

~dt 



(10a) 



Here, t = (l/r)(dH/d<j>) — [— sin^>, cos</>, 8 x cos(j) 
+9 y sm(j)] is the unit tangent vector to the ring (to 
first order in 9), s is the distance along the ring cir- 
cumference measured from, say, the ascending node 
where h(<fi,t) = 0, and 



dR 

~dt 



d6 x . d6 v 
0, 0, r— — sm0 — r—^-cosq 

(XL (XL 



From equation (9) we have 



6L, = \irU ne y • 



SL y = Mr 2 



2 dt 



The torque on the ring is 

c 2ix 



-s 

JO 



rd4> R x F 



or 



T x = 



f 

Jo 

f 

Jo 



(rsin(j>F z — hF y ) 



(rcos(f)F z - hF x ) 



(106) 

(11a) 
(116) 

(12a) 

(126) 
(12c) 



The ring is assumed to exert no torque on itself. Con- 
sequently, the force F (per unit circumference) which 
contributes to T in equation (12a) is due to the inner 
disk and the halo. Thus we have 



-{M/2irr)K 2 z z , 



where 



2,2 

K idz + K hz > 



(13a) 

is the vertical epicyclic frequency excluding the contri- 
bution due to the outer disk [see equation (3)]. Also, 



F x = - (M/2nr)ti rcosq 
where 

n 2 = Ql 



F v = -(M/2nr)n rsincj) 



■ni 



which excludes the outer disk contribution to f2 2 [see 
equation (2)]. Thus we find 

T x = -±Mr 2 (K 2 z -n 2 )9 x , (14a) 



Ty = - l -Mr 2 {K 2 z -tf)6 y . (146) 
The equations of motion d(6Li)/dt = T are 

= -l(n 2 z -n 2 )e x , (15a) 
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- 

dt 


V dt 


d 


( My 


dt 


V dt 



9y , (156) 



after multiplying by two and letting / = Mr 2 denote 
the moment of inertia of the ring. The terms of the 
left-hand side oc are due to the Coriolis force. 

It is useful to introduce 



e(t) = e x {t) - i6 y {t) - 9(t)cx P \-itp {t)j 



(16) 



This representation is well-known from treatments of 
spin precession in quantum mechanics. We can then 
combine equations (15a) and (15b) to obtain 



-A 2 9 



where 



A 



if 



(17) 



(18) 



Note that we can have A 2 < for a prolate halo. 
Solutions of equation (17) can be taken as 6 = 
Cexp (— iu>t) with C a complex constant. From equa- 
tion (17) we have the dispersion relation (w — O) 2 = 
fl 2 + A 2 which gives 



n ± v 7 ^ 2 + a 2 . 



(19) 



In all cases considered here, ft 2 + A 2 > 0. For £1 > 
0, the plus sign corresponds to the fast mode, lo 2 > 
0, with prograde precession [ <p{t) increasing ], and 
the minus sign to the slow mode, u)\ , with retrograde 
precession [ ip(t) decreasing ] (HT). For |A 2 | <C Vl 2 
and f2 > 0, the slow mode has u>\ ~ —A 2 / (20). These 
two modes are the analogues of the normal modes of 
vibration of a non-rotating mechanical system. That 
there are two modes rather than one for each ring is 
due to the Doppler splitting from the rotation (f2). 
A general solution of equation (17) is 

G(t) = Ciexp(-iwif) + C 2 exp(-iw2t) , (20) 

where u>i and loi are the two roots given in equa- 
tion (19), and C\ and C2 are complex constants. The 
four real quantities determine the initial values ^(0), 
y (0), d0 x (0)/dt, and d0 y (0)/dt. If only the fast or 
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slow mode is excited (C\ or C2 = 0), then 9 = 
Const, and ip = u a t+ Const.' with a = 1,2. A 
simple way to excite only the slow or fast mode is 
to take 9 X {0) = C a , 9 y (0) = 0, d6 x (Q)/dt = 0, and 
d9 y (0)/dt = LL> a C a . Figure 4 illustrates the two limit- 
ing cases of Ci /C2 <C 1 and C1/C2 3> 1 . 

We can recast equations (15) in a form equivalent 
to that given by HT by multiplying equation (15a) 
by sin</> and equation (15b) by coscj) and subtracting 
the two. Noting that dh/dcj) = r9 x cos<p + r9 y sin(j) and 
d 2 h/d(j) 2 = —h, we obtain 



d_ 

dt 



n 



d 



(*2 



n 2 



n 2 )h. 



Note that h(r, <j), t) is non-zero only for r equal to the 
ring radius. From equation (2) we have 

n 2 -Tf = . 

r dr 

For a toroidal ring of minor radius Ar <C r, this dif- 
ference is of the order of GMln(r / Ar)/r 3 . In the 
absence of a halo, the right hand side of the equation 
for h is 

-( K 2 dz + n 2 -n 2 )h, 

which is the same as the HT expression 



v S(rQ [h (r,cj>,t)-h(r',<j>',t)] 



• — r'l 3 



The K 2 dz term comes from the integration of over r' < 
r while the fi 2 — fl term is from the integration over 
r' i=a r. 

Multiplying equation (15a) by d9 x /dt and (15b) by 
d9 y /dt and adding the two gives d£ r i ng /dt — 0, where 
the energy of the tilted ring is 

£rin g = \l (el + 91) 2 + \m 2 {9 2 x + 9 2 y ) = Const. , 

(21) 

where a dot denotes a time derivative. This expres- 
sion agrees with the result of HT (their Appendix 
B). In the absence of a halo or with an oblate halo 
(k 2 1z > fi^), we have A 2 > 0, so that both terms in 
Sring are non-negative which implies stability of the 
ring tilting. In the case of a prolate halo with A 2 < 
instability is possible. 

In order to understand the stability of the ring tilt- 
ing, it is useful to consider the influence of a small 
Newtonian drag or friction on the ring motion, a force 



F-^ = — (M/2Trr)(3dIl/dt\ s per unit circumference, 
with (i > but 1 < fl 2 . This drag could result 
from dynamical friction of the ring with the dark halo 
matter (Nelson and Tremaine 1995). 

The full calculation of (3 due to dynamical fric- 
tion on a precessing ring is complicated (Nelson and 
Tremaine 1995), but an estimate based on treating the 
ring as two point masses M/2 at the points max|/i| 
(Weinberg 1985) gives 

13/9, - 2fGM£nA/(3 v ^v 2 r) 



0.035feA 



M \ / 200km/s \ 2 / I0kpc \ 



(22) 

where £nA = 0(3) is the Coulomb logarithm, v Q is 
the circular velocity, and / = r 2 (r 2 + 3r 2 )/(r 2 + r 2 ) 2 
accounts for the density profile of the halo, with r the 
halo core radius. Later, in §3, where we represent the 
outer disk as N tilted interacting rings, the pertinent 
estimate for (3j for a single ring is this formula with 
M the total mass of the outer disk rather than the 
mass of a single ring (Weinberg 1985). 

The drag torque due to F^ is T/ = -\Mr 2 (3 
(d9 x /dt) and T/ = -\Mr 2 {3{d9 y /dt). Solution of 
equation (17) including T-^ gives 



= n ± Vn 2 + a 2 - — (i± 



n 



Vn 2 + a 2 



• (23) 



Thus, both fast and slow precession modes with 
A 2 > are stable in the presence of friction, whereas 
the slow precession mode with A 2 < is unstable 
[ LUi = Im(oj) > ]. For the slow precession mode, 
the kinetic energy term cx 9 2 y in Sring is smaller than 
the magnitude of the potential energy term cx 9 2 
for I A 2 1 < ft 2 . In the presence of the friction force, 



d£. 



,/dt = -I[3{9 2 x + 91). Thus, for the slow pre- 



""-rings 1 1 v y/' ^w.. 

cession mode, £ r ing oc A 2 is negative if A 2 < 0, and 
the tilt angle grows. The instability is therefore a 
negative energy dissipative instability. An analogous 
instability was predicted and observed for the pre- 
cession of laboratory collisionlcss rclativistic electron 
rings (Furth 1965; Beal et al. 1969). The instabil- 
ity also occurs in the presence of a non-Newtonian 
(nonlinear) friction force on the ring. 

A second quadratic integral of the motion of equa- 
tions (15), analogous to angular momemtum (see fol- 
lowing), can be obtained by multiplying equation 
(15a) by ~9 y and equation (15b) by 9 X and adding. 
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This gives dV r i ng /dt = 0, where 

Vring = I (0 y 6 x ~ 9 x 9yj + (9% + 9y) , 

= I9 2 (n-ip) = Const. (24) 

An analogous constant of the motion exists for cases 
of two or more rings (see §§2.3, 3, and 4). This con- 
stant of the motion is a new result of the present 
work. For the case of instability, with u> complex, 
d<p/dt = ut r = Re(uj). With the above-mentioned 
friction force included, we have dVring/dt — [3I9 2 ut r , 
which is compatible with equation (23) which shows 
instability for A 2 < and gives < uj r < £1 and 

V rings ^ 0* 

The second integral (24) can be understood by not- 
ing that equations (15) follow from the Lagrangian 

£ = \i(fi + 8i) +m (0y0 x -6j v ) 

-I/A a («£ + «;) . (25a) 

A canonical transformation to the variables 9 and <p 
gives 

C = I J (e 2 + 0V) - IQ9 2 ip - \lA 2 9 2 . (256) 

Thus the canonical momentum P v = dC/dip = I9 2 
(ip—Cl) = Const. = —Vrings because d£/dip = 0. On 
the other hand, Pg — 19 depends on time. The Hamil- 
tonian H(0,<p, P e , P v ) = 9(dC/d9) + ip(dC/dtp) - C 
is 

and H — £ r ing = Const, in that £ has no explicit 
time dependence. 

The dependent terms of H can evidently be 
viewed as an effective potential for the 9 motion. In 
general (for P v ^ 0), 9 nutates between a minimum 
and a maximum value. The angular frequency of the 
nutation is 2(fi 2 + A 2 ) 2 so that the nutation period is 
Pnut = CP/2)(1 + A 2 /n 2 )-?, where P = 2tt/0 is the 
period of the orbit. This can also be seen by noting 
that with 0(f) given in general by equation (20), we 
have |0(f)| = |Ci+C2exp[— i(u)2— ^i)f]| andw2-^i = 
2(fl 2 + A 2 )i. For |A 2 |/f7 2 = e < 1, the period of the 
slow mode Pi — 2w/\ui\ w (2/e)P S> Pnut, while 
the fast mode period P 2 = 2tt/uj 2 ~ (-P/2)(l - e/4) 



is slightly longer than the nutation period P nu t ~ 
(P/2)(l — e/2). The nutation of the ring is clearly 
evident in Figure 4. For \C 1 /C 2 \ < 1 (Figure 4a), 
0(f) nutates many ( w 4/e ) times in the period Pi of 
motion of about the origin. On the other hand, for 
\Ci/C 2 \ > 1 (Figure 4b), 0(f) has an elliptical path 
with the azimuth of say the maximum of |0(f)| pro- 
cessing slowly in the clockwise direction with a period 
« (2/e)P. 

A counter-rotating ring with — > — O behaves in 
the same way as a co-rotating ring. For A 2 > 0, the 
ring slow mode precession is retrograde relative to the 
ring particle motion but of course prograde relative to 
our coordinate system. 

2.2. Interpretation of Ring Observations 

Warps of spiral galaxies are deduced from measure- 
ments of 21 cm line neutral hydrogen emission. Here, 
we consider the spectral signature of a single ring in 
either the slow precession mode or the fast precession 
mode. The observer is considered to be in the (a;, z) 
plane at an angle to the z— axis. The velocity of 
the ring matter is given by equation (10). Thus, the 
velocity of the HI in the direction of the observer is 

v\\{4>) / {r£l) = — sin(^>)sin(0) 

+ 9 (l - ^) cos(V0cos(^ - <p) ■ (26) 

The spectrum of the ring is iS(um) oc \dv\\/ def)^ 1 

oc (1 — v 2 /v 2 max )~^ ■ The ring motion dH/&t\ 8 en- 
ters through the term involving u> a . For a ring in the 
slow precession mode with |wi/f2| <C 1, the ring mo- 
tion has only a small affect on V|| (</>). The situation 
is very different for a ring in the fast precession mode 
where u 2 /£l > 2 if A 2 > 0. In this case the sign of the 
term in vu (</>) proportional to 9 is reversed which cor- 
responds to a change of the line-of-nodes angle ip by 
±180°. Observations are interpreted assuming neg- 
ligible ring motion, |w a /fi| <C 1, and this is wrong 
if the fast precession mode(s) is excited. The fast 
precession modes of a system with many rings may 
disappear from view over a long enough time due to 
phase-mixing (see §3.1 and §3.2). 

2.3. Two Tilted Rings 

Consider now the case of two tilted plane circular 
rings of mass Mj and angular rotation rate Slj at radii 
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Vj with j = 1,2. We have 



SLj x — — — 7- h Ij£lj6 



2 -3 df T^-j-jy , 
1 d& 



(27a) 



where 6*.^ and 6*j y are the tilt angles and 7j = Mj-r 2 
the moments of inertia for the two rings. 

To obtain the torque on, say, ring 1, we exclude as 
before the force due to this ring on itself. Thus the 
horizontal force on ring 1 is 

(Fi)*,„ - -Mi (n 2 d + n 2 + 1^®) ( Xl , yi ) , 

(28) 

where the <&r 2 is the gravitational potential due to 
ring 2, and the 1-subscript on the parenthesis indi- 
cates evaluation at r = r\. Note that 



7/2 



1 d 



GM 2 



r dr ) 1 n dri J Q 27r |ri — r 2 | 

= GM 2 I 12 - GM 2 r 2 Ji2/ri , (29) 



where 



I12 = 



Jl2 = 



1 



2n [r\ +r\- 2rir 2 cos(f>y 
2n d<h cosd) 



(30a) 
(306) 



h 2lT ( r \ +r\- 2r 1 r 2 cos<p) 2 
There is also the vertical force on ring 1, 

F lz = - (Mi^Trn) (n 2 dz + K L)x fti + Fl 2 , (31a 
where 

r27r #' GM 1 M 2 [h 1 ^,t)-h 2 (4>',t)} 



r n df_ GM 1 M 2 [1 
Jo 27r [r 2 +r 2 - 



2rir 2 cos((j) — </>')]■ 



(316) 

is the vertical force on ring 1 due to ring 2. This 
integral can be simplified by noting that 

h 2 (4>\ t) = r 2 9 2 sm(<l) — ip 2 + <j> — <f>) 

= h 2 cos(4> — 4>) — r 2 9 2 cos((j) — ip 2 )sm(<j) — 4>) . (31c) 

The sin(0 — <fi') term in equation (31c) does not con- 
tribute to the integral (31b). Thus we have 



/ n# Ri x (zff) = 
Jo 



I 1 

--GM x M 2 r\l X2 d Xx -k + -GM x M 2 r x r 2 J X2 e 2x k 

-iGM 1 M 2 r 2 / 12 ^y+iGM 1 M 2 r 1 r 2 Ji2^ !/ y • (32) 
Combining terms in equations (12b) and (12c) gives 

T lx = -^M ir jAle lx - l -C 12 (0 lx - 9 2x ) , (33a) 

Tiy = -^MrfAfay - \c l2 (6 ly - 6 2y ) , (336) 
where 



C12 = GM x M 2 r x r 2 J X2 



(34) 



measures the strength of the gravitational interaction 
between the two rings, and A 2 is defined in equation 
(18). The two terms involving the integral I\ 2 cancel. 
Note that C\ 2 is symmetric between the two rings and 
that it is non-negative because J\ 2 > [ see equation 
(45e) and Figure 5 ]. 

The torque on ring 2 is found in the same way to 

be 

T 2x = -^M 2 r 2 2 A 2 2 e 2x - X -C 12 (8 2x - 9 lx ) , (35a) 



T, 



2y = - l -M 2 r 2 A 2 6 2y - l -C 12 (9 2y - 6 ly ) . (356) 

The equations of motion are 

d{SL{)/dt = Ti , d(SL 2 )/dt = T 2 . (36) 

For the limit where the only torques are those due to 
the two rings, A 2 = and A 2 . = 0, we find 



j f {5U + SL 2 ) = , 



(37) 



which is a necessary result. 

We can follow the steps leading to equation (22) to 
obtain d£ r i ngs /dt = for the energy £ r i ngs of the two 
tilted rings, where 

1 2 1 2 

Srtngs = ^ ^ Ij (0% + 2 y ) + - ^ J J A j (tfx + °jy) 



1 2 1 2 

where I\ = M\r\ and I 2 = M 2 r 2 are the moments of 
inertia of the two rings. This energy is non-negative 
and the rings are stable if both A 2 and A 2 , are non- 
negative. Also, in this case the rings are stable in 
the presence of dissipative forces such as the friction 
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force of Subsection 2.1. Thus a necessary condition 
for instability with or without dissipation is that A 2 
and/or A 2 be negative. 

A second quadratic integral of the motion can be 
obtained by following the steps leading to equation 
(24). We find dV rtngs /dt = 0, where 



/l0?(fil-¥>l) + / 2 02(ft2-¥>2) • (39) 



For the case of instability uj — uj r + iui, uoi > 0, we 
have dipj/dt = uo r . Thus, a further necessary condi- 
tion for instability in the absence of dissipation (in 
addition to A\ and/or A 2 , being negative) is that 
ill — oj r and £1 2 — <jj r have different signs; that is, 
Q, 2 < uj r < ^1 for 0,2 < fii- For the slow precession 
modes with \oj r \ <C ilj, this corresponds to counter- 
rotating rings, Q1/CI2 < 0. 

The Lagrangian for the two ring system is 

3=1 

\c l2 [e\ + e\ - 20 1 e 2 coB{<p 1 - . (40) 

The corresponding Hamiltonian H — ^{OjdC/dOj + 
ipjdC/dipj) — C is equal in value to £ r ings of equation 
(38). The second constant of the motion, V rings = 
—P v \ — P V 2, where P V j = d£/d(pj. Note that 
dP vl /dt = dC/dif! = 
dP v2 /dt = d£/dtp 2 
d{P v i + P v2 )/dt = 0. 

As in the case of a single ring, it is useful to intro- 
duce 

Q 1 (t) = e lx -ie ly = e 1 (t)exp[-ii Pl (t)] , (41a) 

e 2 (t) = e 2x - ie 2v = e 2 (t)cx P [-^(t)] . (416) 

The equations of motion (36) can then be written as 

C12 



|#i# 2 sin(<£>i — 1P2) and that 
2#i#2sin(</?i — ip 2 ) so that 



(6i - e 2 ) 



(42a) 



-A 2 2 -^(0 2 -0l) 



/ 2 



(426) 

after dividing through by Ij/2. 

For ©i and 2 oc exp(— iut), we get the dispersion 
relation 



(w - Oi) - # 2 (w - n 2 ) - 



°12 
hh 



(43) 



where 



1 3 



with j = 1,2. Equation (43) can readily be solved for 
the four roots. Two of these roots are high-frequency, 
fast precession modes with u> > 2f2j , assuming Oj > 
and A 2 > 0. The other two roots are the slow preces- 
sion modes with much smaller frequencies, typically, 
lu 2 « <>;. 

In order to get simple analytic results, we first con- 
sider the low- frequency modes lo 2 -C £1 2 , assuming 
|A 2 | < tt 2 and Ci 2 //j < ft 2 . We can then neglect the 
second time derivatives in equations (42) and the uo 2 
terms in equation (43), and this leads to a quadratic 
dispersion equation for u>. The roots of this equation 
are 

u = 2 t^P 1 + w p 2 _ w ff 12 ^ + ^ 12 )1 



± 



UJpl + LOgl2 (1 - 4 2 )] 2 + 4Wg 12 ^l 2 | 2 . 

(44) 



Here;. 



2Qi 



(45a) 



for j = 1,2, are the slow precession frequencies of the 
two rings in the absence of gravitational interaction 
between them (C12 = 0); 

u g i2 = ^ and £12 = ^- , (456) 
IL\ L2 

where lo q \2 is a frequency which measures the strength 
of the gravitational interaction between the two rings, 
and where Lj = Ijflj are the angular momenta of the 
two rings. We have 



0g\2 



2.8 x 10" 



M 2 



10 10 M, 



20 kpc\ /200 km/s 
r 2 ) \ Oin 



i 2 > 



where 



12 



2tt 



5/ ^ 



COSfl 



2?r [1 + S 2 - 25cos0] 3 



(45c) 



(45e) 



with <5 = n/r 2 < 1. For 5 < 1, J{ 2 -» §6 2 . For 
1 - 5 < 1, J( 2 -» 1/[tt(1 - <5) 2 ]- We assume that 
1 — 6 > Az/r, where Az is the disk thickness. The 
dependence of J{ 2 on 5 is shown in Figure 5. 
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A simple limit of equation (44) is that where A 2 = 
= A 2 . The two roots are then w = and cj = 
-w fl i 2 (l + M = -\C X2 L X L 2 I{L X + L 2 ). The zero 
frequency mode corresponds to both rings tilted by 
the same angle, 9i = 6 2 , which is the rigid tilt mode 
(HT). The other mode has 9i = -(L 2 /ii)6 2 ; that is, 
the rings are tilted in opposite directions for L\jL 2 > 
0. 

We first consider co-rotating rings where the four w 
roots of equation (43) are real in the absence of dissi- 
pation. Note that the low frequency roots of equation 
(44) are real for €12 > 0. A general solution for the 
motion of the two rings is then 



@j = ^ Cj a exp(—ioj a t) 



(46) 



where {uj a } are the four mentioned frequencies, and 
{Cj a } are 8 complex constants. However, only 8 real 
quantities are needed to specify general initial con- 
ditions because the {C\ a } are related to the {C 2a } 
by equations (42). If only a single mode of the 
system is excited, then 0i = Ci ( gexp(— iujpt) and 
02 = C2f3cxp(—iu)pt). Equations (42) then imply 
that 62/61 is real. Because Qj = 9jexp(—iipj), we 
conclude that 62/61 oc exp[— i(<p2 — fi)] is real so 
that if2 — <fii is either or ir. That is, the rings 
are either tilted in the same direction and precess to- 
gether or they are tilted in opposite directions and 
also precess together. From equations (42), 62/61 = 
h [Kl - (up - r!i) 2 ] /C12. 

Figure 6a shows the u> g i2 dependence of the fre- 
quencies of the two slow-precession modes obtained 
from equations (43) and (44) for the case where both 
rings are co-rotating and the halo is oblate. The 
'slower' of the two modes has 62/61 > 0, while the 
other has 62/61 < as shown in Figure 6b. The 
larger gravitational torque in the second case accounts 
for the faster retrograde precession. As uj g i2 increases, 
the frequency of the slower mode approaches 



L\ + L2 L\ + L2 



(47) 



and 62/61 — > 1 for uj gl2 > \uJ p2 - w p i|/(l +^12), 
where the ui p j are defined in equation (45a). This re- 
sult can be derived from equation (43) or (44). For 
the case of N rings and sufficiently large tv g jk, we 
find in general that the lowest frequency approaches 
the angular momentum weighted average of the sin- 
gle ring precession frequencies u) p j [equation (45a)]. 



In contrast with the dependence of u>i, L02 decreases 
monotonically as uj g i2 increases. 

From the approximate equation (44) for the slow 
precession modes, we can distinguish weak and strong- 
coupling limits of two co-rotating rings. The coupling 
strength is measured by the dimensionless parameter 



?12 = 



+ tl2) 



\uj p 2 — (dpi I 



GJ12 (Li + L 2 ) 



2(f2in)(f2 2 r 2 )|a;p 2 - u p i\ ^ 

which is symmetrical in the ring indices (unlike w s i 2 ). 
For weak coupling, £i 2 -C 1, the two rings are affected 
little by their gravitational interaction, and the two 
slow mode frequencies arc uu\ ~ iv P 2 with 62/61 3> 1 
and LO2 ~ Lo p i with I61/62I 1. In the opposite 
limit of strong coupling, £i 2 3> 1, the gravitational 
interaction of the two rings is important, and u)\ is 
given approximately by equation (47) with 62/61 ~ 
1, and \uj 2 I > |wi| with 6 2 /6i < 0. 

For nearby rings, (Ar) 2 = (r 2 — ri) 2 <C r\ , of equal 
mass Mi = M2 = 2-irrAr S (with E the surface mass 
density of the disk at r), we have Ji 2 w l/[7rr(Ar) 2 ], 
\uj P 2 — u) p i\ w \du p /dr\Ar, and thus 



4 



2GS 



12 



ft(Ar) 2 \duj p /dr\ 



(49) 



Evidently, there is strong (weak) coupling for Ar <C 
Ar c (Ar ^> Ar c ), where Ar c corresponds to ^i 2 = 1 
so that 

2GE 



Ar c = 



(50a) 



Jl\duj p /dr\ 

Assuming a flat rotation curve (fir = const.) and 
(A/fl) 2 — const., we obtain 



Ar c w 13.2 kpc 



200 km/s 
Or 



20 kpc 



0.2 



-V 



|A 2 |/r» 2 ; VioMo/pc 2 / (506) 

Rings closer together than Ar c will be strongly cou- 
pled. 

Figure 7 shows sample orbits of 61 and 6 2 for two 
co-rotating rings in an oblate halo for initial condi- 
tions 6i/6 2 = 1/2 and G»i = = 6 2 . The top panel 
shows a case of weak coupling (small £12) where the 
line-of-nodes <pj of the two rings regress essentially 



9 



independently with the result that ip2 — fi increases 
linearly with time. The bottom panel is for a case 
of strong coupling (large £12) where the line-of- nodes 
are phase-locked in the sense that the time average of 
<P2 — <Pi is zero. Figure 8 shows the azimuthal angle 
difference of the line-of-nodes for the two cases of Fig- 
ure 7. In the case of weak coupling, a low frequency 
initial perturbation excites both low frequency modes 
which have similar frequencies (oj 2 ~ 2o>i), and this 
gives the unlocked behavior. For strong coupling, the 
initial perturbation excites mainly the u>i mode be- 
cause \oj2\ ^> |wi|, and this gives the locked behav- 
ior. For the parameters of figures 4 and 5, note that 
£12 ~ 60a> s i2/f2i. 

Figure 9a shows the uj g i 2 dependence of the fre- 
quencies of the two slow-precession modes obtained 
from equation (43) for the case where both rings are 
co-rotating and the halo is prolate. (We have adopted 
the convention of ordering the mode frequencies by 
their magnitudes for ui g i2 — > 0.) Figure 9b shows 
that the lo 2 mode has 62/61 > 0. As in the case of 
an oblate halo, for large u g \ 2 or £12, u>2 approaches 
the limit given by equation (47) and 62/61 — > 1. 

As in the case of a single ring in a prolate halo 
potential [see equation (23)], small friction torques, 
with friction coefficients /3i and (3 2 , on two co-rotating 
rings in a prolate halo lead to a negative energy dissi- 
pative instability. We find that the slow modes with 
prograde precession, u> a > 0, are unstable in the pres- 
ence of dissipation. That is, the lo\ mode in Figure 9 is 
unstable (stable) for uj gl2 < (>) 0.02, whereas the w 2 
mode is unstable for all u> g i2- The two high frequency 
modes are prograde but are damped by friction for all 
uj g i2- As in the case of an oblate halo (Figure 8), we 
observe phase-locking for sufficiently large values of 
£12. 

Frictional torques between two co-rotating rings in 
a prolate halo also give rise to a negative energy dis- 
sipative instability. Such torques could arise from the 
viscous interaction between adjacent gaseous rings. 
These frictional torques do not change the total an- 
gular momentum of the two rings. Therefore, these 
frictional torques, if Newtonian (linear) in nature, are 
necessarily of the form that they contribute terms 




on the right hand sides of equations (42a) and (42b), 
respectively, where (3\2 is the friction coefficient (with 
units of angular frequency). For a gas with kinematic 
shear viscosity v (with units of cm 2 / s) the terms (51) 
arise from the momentum flux density or stress T rz 
due to the different vertical velocities of the two rings 
as a function of <j>. This gives (3\2 = ^/(Ar) 2 , where 
Ar = T2 — r\. The relevant viscosity is that due to 
turbulence in the gas v = v t owing to the smallness 
of the microscopic viscosity. An estimate of v t can be 
made following the proposal of Shakura (1973) and 
Shakura and Sunyaev (1973) that v t — ac s Az, where 
c s is the sound speed, Az is the half-thickness of the 
disk, and a is a dimensionless constant thought to be 
in the range 10~ 3 to 1. In Sections 3 and 4 we discuss 
further the frictional torques between adjacent rings. 

2.4. Two Tilted, Counter-Rotating Rings 

Here, we consider the case of two tilted rings which 
are rotating in opposite directions. This situation is 
pertinent to observed counter-rotating galaxies (see, 
for example, Jore, Broeils, and Haynes 1996). From 
equation (44), instability is possible (in the absence 
of dissipation) only for £12 < in agreement with our 
discussion of equation (39). Thus, consider ring 2 to 
be counter-rotating (f2 2 < and L 2 < 0), while ring 
1 is co-rotating. The region of instability is bounded 
by the two curves 

A 2 A 2 

UJp2 ' UJpl -Wh\ + Wh\ 

= (-l±2|*i 2 |'-|*l 2 |)w fl 12, (52) 

on which the square root in equation (44) is zero. 
The curves are shown in Figure 10. For parameters 
in the region between the two curves there is insta- 
bility, which is seen to occur only for A 2 and/or A 2 
negative. For fixed values of u> p \, lu P 2, and £12, in- 
stability occurs in general (\£i2\ 7^ 1) for w 5 i 2 larger 
than a critical value but smaller than a second larger 
critical value. The maximum growth rate is for 

(LU p2 - LU p l)/uj g l2 = -1 - |^12 1 (SO that LU p2 < LOpl), 

and it is max(wi) = L0 g \2\l\2\^ ■ This is a dynamical 
instability in that it does not depend on dissipation 
as in the case of a single ring with A 2 < 0. For con- 
ditions of maximum growth, 61/62 = — %l\lv^, so 
that 9i/9 2 = \L 2 /Li\i and ipi — ip 2 = 7r/2. We can 
thus view the co-rotating ring as 'torquing up' the 
more slowly precessing counter-rotating ring by the 
gravitational interaction or vica- versa. 
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Figure 11a shows the Lo g i2 dependence of the fre- 
quencies of the two slow-precession modes obtained 
from equation (43) for the case where ring 2 is 
counter-rotating (O2 < 0) and ring 1 is co-rotating 
(fli > 0) and the halo is oblate. Figure lib shows 
that the mode w\ with prograde precession involves 
mainly motion of ring 2, whereas the 0J2 mode with 
retrograde precession involves mainly ring 1. For the 
conditions studied (w 9 i2 < O.lili) , we find no phase- 
locking. 

2.5. Three Tilted Rings 

Extension of the results of Subsection 2.2 to the 
case of three tilted rings is straightforward. Figure 
12a shows the dependence of the three slow precession 
frequencies on the strength of the gravitational in- 
teraction measured by w g i2 with ui g 23 / oj g \2 = Const. 
For increasing w g i2, the slowest mode (wi) slowly ap- 
proaches the limiting frequency given by the gener- 
alization of equation (47) to three rings. The other 
two frequencies decrease monotonically with increas- 
ing w g i2- Figure 12b shows the geometrical nature 
of the three slow precession modes. Figure 13 shows 
sample orbits &j(t) for three rings for cases of weak 
and strong coupling. Figure 14 shows different cases 
of phase-locking for a three ring system. 

3. N Tilted Corotating Rings 

The generalization of equation (40) gives the La- 
grangian for an N ring system as 

1 N 



N 



~ \ Y. C Jk W 3 +0l- 20 7 -0 fe cos(<^ - y k )] . (53a) 
j<k 

The Hamiltonian, 

dC dC 

^ d0 k otpk 



is 



i=i L 
1 N 



f + +<»,)' 



jk 



fj +6l- 26» j 6» fe cos(^ - ip k ) 



(536) 



Because the {C^} are non-negative, all of the terms 
in H are non-negative if A| > which is the case for 
oblate halo. Thus the N— ring system with A 2 > is 
stable in the presence of dissipative forces. The total 
canonical angular momentum is 



r>tot 



N N 



^Itftej-ilj). (53c) 



In the absence of dissipation £ r ings = "H and V r ings = 
— P^ ot are constants of the motion. 

As in the case of one or two rings, it is useful to let 

9,(t) - 6 jx (t) - iO jv (t) = ^(t)exp[-*^(t)] . (54) 

The equations of motion can then be written as 

N 

-ApjQj-^CjkiQj-Qk), (55) 

with j = 1..N, where the Cjk generalize C12 of equa- 
tion (34). 

As discussed in §2.1, a Newtonian drag force on 
the ring motion due to dynamical friction with the 
halo matter can be taken into account by including 
the term 

,/<-> ■ 

(56) 



-8lA 

Pj 3 dt 



on the right hand side of equation (55). The influ- 
ence of the relative friction discussed in §2.3 can be 
accounted for by including the terms 

-^• + iJ J - + i(e J -e j+ i)-^._i/ J -_ i (e J -e j _i) (57) 



on the right hand side of equation (55). Here, Ij±i = 
(/,-Jjii)* and 0'. ±1 = f3j j±1 



^ t /(Ar) 2 ] J+ i with 



0\ = and 0' 



N+4 



0. With the dissipation terms 

(56) and (57) included one can readily show that 
dSrings/dt < (see also §4). We have not found an 
analogous result for dV r ings/dt in the presence of dis- 
sipation. 

We have developed and tested codes to solve equa- 
tions (55) (including the terms (56) and (57)) for 
N — 7, 13, 25, and 49 rings. For the results pre- 
sented here, the rings are taken to be uniformly 
spaced in r with ring one at n = 10 kpc and ring 
N at r N = 40 kpc so that rj = 10 + (j - l)Ar, 
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with Ar = 30/ (N — 1) in units of kpc. The ring 
masses Mj = 27rr 7 -Ar£(rj) and moments of inertia 
Ij = MjTj are calculated with £(rj) the sum of the 
surface density of the inner exponential disk and that 
of the neutral hydrogen given in §2. The A| are calcu- 
lated using equations (4), (6), and (18). The coupling 
coefficients {Cjk} are evaluated using equations (34) 
and (30b) and stored. Equations (55) are solved as 
four first order equations for 6 x j, 6 x j, 9 y j, and 9 y j. 
At the same time, an additional equation, 



o2 I al 
xj "r" yj 



(58) 



is also solved to give ipj (t) , which is the line-of- nodes 
angle relative to the x— axis. 

Comparisons of the temporal responses of 25 and 
49 ring sytems for different initial conditions show 
generally good agreement of the warps 8j for the 
time intervals studied. This indicates that N = 49 
rings gives a valid representation of a continuous disk. 
However, in contrast with a continuous disk (HT), an 
N— ring disk has 2N discrete modes with separate 
frequencies - N 'low-frequency' modes and N 'high- 
frequency' modes. Figure 15 shows the power spec- 
trum of 6 X 2b{t) for a 49 ring disk-halo system which 
is initially given a random perturbation with all ra- 
dial wavelengths excited and no damping. The fre- 
quencies cj are normalized by the angular velocity of 
the inner ring Ql Note that Pi = 27r/Oi w 0.308 
Gyr for our reference values. The frequency differ- 
ences between adjacent modes is Su) ~ Qi/N. For 
treatment of initial value problems of disk warping 
for time intervals < t max , frequency differences 6u> 
less than l/t max are irrelevant. It is sufficient to 
have Slo ~ l/t max which corresponds to N ~ Q\t max 
= 27r(t maa; /0.3Gyr) rings. This suggests that a larger 
number of rings may be needed for our cases with 
tmax > 8 Gyr (Figures 18 and 19). This will be in- 
vestigated in a future work. 

The lowest frequency mode of the power spectrum, 
uj\ in Figure 15 (see also §2.3 and §2.5), remains dis- 
crete (that is, isolated) as N — > oo because our disk 
model has a sharp outer edge so that 1/S(r) is inte- 
grable (HT; Sparke and Casertano 1988). Power spec- 
tra obtained for cases with Cjk — * "(Cjk show that u)\ 
increases with 7 much more slowly than the higher 
frequency modes which is consistent with the behav- 
ior observed in Figures 6 and 12. For the conditions of 
Figure 15, the magnitude uj\ is considerably smaller 
than that predicted for the limit of strong self-gravity 



between the rings where u>i — > "}2 LjU p j / "^2 Lj [equa- 
tion(47)]. This limiting frequency is the same as the 
frequency u t of Sparke and Casertano (1988). In any 
case, the frequency ui\ is so small (the corresponding 
period is > 8 Gyr) that it is irrelevant to the initial 
value problems considered below. 

In the following four subsections we consider dif- 
ferent possible origins of warps in an otherwise flat 
galaxy. Subsection 3.1 discusses warp excitation by a 
passing compact satellite; §3.2 treats warp excitation 
by a compact sinking satellite; §3.3 treats the case of 
a tilted halo potential. Subsection 3.4 considers the 
tilt evolution for the case where the initial plane of 
the disk material is tilted. 

3.1. Warp Excitation by a Passing Satellite 

Consider the excitation of a warp in a galaxy due 
to the passage of a satellite of mass M s much less 
than the galaxy's mass. The satellite's orbit r s (t) 
can easily be calculated exactly for M s — > 0, neglect- 
ing the back reaction of the perturbed galaxy on the 
satellite. However, we assume that even at closest 
approach the satellite is far from the center of the 
galaxy and therefore calculate the orbit in the halo 
potential, equation (6a), with the ellipticity of the 
halo neglected for simplicity (q = 1). At the closest 
approach at min(|r s |) = r so , the satellite's speed is 
v so , and the angle between the satellite and the (x, y) 
plane is 6 so . Also, at closest approach, the satellite is 
taken to be at (r s ) y = and (r s ) z > 0. We consider 
both prograde and retrograde satellite passages. 

In the presence of a satellite of finite mass M s , the 
description of the galaxy disk in terms of tilted cir- 
cular rings breaks down. The centers of the rings are 
shifted from the origin, and the rings become non- 
circular. For this reason we consider the response of 
the galaxy to a 'symmetrized' perturbation of a satel- 
lite obtained by replacing the actual satellite by two 
satellites, each of equal mass M s /2, with one satellite 
in the orbit r s (t) described above and the other in the 
orbit — r s (i). With this prescription we then calculate 
the torque of the satellites on each ring as a function 
of time, 



' W 2 3 J 2ir 



+ ( r s 



d * t 

- — $ s {Tj,T s 

), 



where Vj(t) — [rjcoscj), r^-sin- 
given by equation (8a), is the position vector to a 



(59) 

<f>, hj(<f>, t)}, with hj(4>, t) 
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point on the jth ring, and j = 1, .., N runs over the 
different rings. Here, $ s = — GM s /[(y } ■ — r s ) 2 + a 2 ] 2 is 
the gravitational potential of the satellite and a is its 
characteristic radius. To obtain the warp response of 
the galaxy we evaluate equations (59) numerically for 
j = 1,..,N and include the complex torques T?(t) = 
Tj x (t) — iTj (t) on the right hand sides of equations 
(55). 

Figure 16a shows the Briggs plot of the warp re- 
sponse 9(<p) of the galaxy resulting from the retro- 
grade passage of a satellite of mass M s = 2 x 1O 1O M . 
Note that the dependence of 6{ip) at the times shown 
is of the form of a leading spiral wave and that it is 
qualitatively similar to the leading spiral waves ob- 
served in many warped galaxies (Briggs 1990). Note 
also that the form of the spiral wave is preserved for 
a long time. However, the amplitude of the warp 
9 is significantly smaller than the warps observed in 
many cases, and this agrees with the conclusions of 
HT. Nevertheless, it is of interest to understand the 
behavior. Figure 16b shows the warp response for 
the retrograde case with and without gravitational 
interactions (Cjk) between the rings. Clearly these 
interactions give a strong phase-locking of the inner 
rings of the disk (roughly, n to r 14 or 10 to 18 kpc) 
which have the same ipj's as a function of time. The 
phase-locking of two and three rings was disscussed 
in §2.2 and §2.4. 

Figure 17a shows the Briggs plot for the warp re- 
sulting from the prograde passage of the satellite. For 
early times (t < 2 Gyr) the form of 0(<p) is a trailing 
spiral wave which is distinctly different from that for 
the retrograde passage at a similar time. The am- 
plitude of the warp is noticably larger than for the 
retrograde case but still smaller than many observed 
warps. For later times (t > 4 Gyr) the trailing wave 
'unwraps,' evolving into a leading spiral wave of form 
similar to that for the retrograde passage for radial 
distances out to about 35 kpc. For larger r, the 0(<p) 
curve has a 'spur' where 9 increase with r while ip re- 
mains roughly constant (see §3.2). This unwrapping 
appears to be due to the phase-mixing of the fast pre- 
cession modes which propagate over the radial extent 
of the disk in times less than 4 Gyr. 

The different behavior of the prograde and retro- 
grade cases for early times (t < 1 Gyr) can be under- 
stood by considering the response of a single ring with 
or without gravitational interactions between rings. 
Figure 17b shows the orbits [9j x (t), 9j y (t)] of ring No. 
25 (radius r 2 5 = 25 kpc) including the gravitational 



interactions for prograde and retrograde encounters. 
Consider the torque on the j th ring Tf = T? x — iT- y 
which enters on the right hand side of equation (55). 
The tilting of the ring during the encounter is rel- 
atively small and can be neglected for this discus- 
sion. Consequently, T? {t) is an even function of t 
if the pericenter is taken as t = 0. Further, T? y (t) 
is the same for prograde and retrograde encounters. 
For the considered geometry of the satellite orbit, 
Tj y (0) < 0, and the magnitude of T? (i) decreases 
monotonically with increasing \t\. On the other hand, 
the £— torque Tj X (t), although smaller in magnitude, 
is an odd function of time, and it has opposite signs 
for prograde and retrograde encounters. Very roughly 
T? x (t) oc d[T? (t)]/dt for the prograde case whereas 
T? x (t) cx -d[T? (t)]/dt for the retrograde case. As a 
function of time, T? traces out a closed clockwise path 
in the complex plane for the prograde case, whereas 
this path is counterclockwise for the retrograde case. 
For the prograde case, the path of T? is in the same 
clockwise sense as the path of the orbit Qj = 9 X j—i9 y j 
for the fast precession mode (sec §2.1), whereas for the 
retrograde case the path of Tf is in the same direction 
of motion as Qj for the slow precession mode. For this 
reason a prograde encounter couples strongly to the 
fast precession mode, whereas a retrograde encounter 
couples mainly to the slow precession mode. 

3.2. Warp Excitation by a Sinking Satellite 

Here, we discuss the behavior of warps excited by a 
slowly sinking compact minor satellite. The satellite 
is assumed to be minor in the respect that its mass is 
much smaller than that of the galaxy plus halo mat- 
ter within say 30 kpc. The sinking of a more massive 
satellite is likely to cause substantial thickening of 
the disk (Walker, Mihos, and Hernquist 1996) which 
is not observed. Initially, the satellite is assumed to 
be in an approximately circular bound orbit in the 
halo gravitational potential at a large radius r so in the 
(a;, z) plane at an ange 9 so above the (x, y) plane. The 
satellite slowly sinks owing to dynamical friction with 
the halo matter [equation (6a) with q = 1] which in 
the simplest description (BT, p. 428) causes the spe- 
cific angular momentum of the satellite l s ~ r s {t)v 
to decrease as 

^± ~ _ .428^f ink , (60a) 
dt r s (t) v ' 
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where £nA is the Coulomb logarithm. The time for 
the satellite to sink from r so to the galaxy's center is 



t sin k « 4.40 Gyr 



1O 1O M 



50 kpc 



200 km/s / \£nA 



(606) 



The disk response is found by integrating equations 
(55) including the numerically evaluated torques (59) 
for r s (t) at each time step for a 'symmetrized' sinking 
analogous to the approach discussed in §3.1. 

Figure 18 shows the Briggs plots of the warp re- 
sponse 8(ip) of the galaxy at different times after the 
sinking (r s = 0) of a satellite of mass M s — 10 10 M Q in 
retrograde and prograde orbits. A number of points 
are observed: 

(1) The warp amplitudes are larger than for the 
passing satellite (of mass M s = 2x 10 10 M Q ) discussed 
in §3.1, but the amplitudes are still smaller than some 
observed warps (Briggs 1990). 

(2) Coiling in the 9((p) curves at early times (t < 2 
Gyr) tends to disappear at later times (t > 4 Gyr). 
We believe but have not proven that this is due to 
the phase-mixing of the fast precession modes which 
propagate over the radial extent of the disk in times 
< 4 Gyr (see also §3.1). 

(3) For both the retrograde and prograde cases 
6(ip) exhibits a leading spiral wave for t > 4 Gyr qual- 
itatively of the form observed for M 83 (see Figure 2c). 
However, the tilt amplitude 9 is significantly smaller 
than that for M 83. 

The dissipative torques due to dynamical fric- 
tion [equation (56)] and that due to relative friction 
[equation (57)] have different consequences for the 
warp evolution. For a dynamical friction coefficient 
= 0.1, the warp amplitudes 9 for the cases 
of Figures 18a and 18b are reduced by factors ^0.5 
while the line-of-nodes angles are roughly the same. 
On the other hand, for relative friction coefficients 
Pj/^i < L short radial wavelengths are damped out 
without appreciably affecting the overall shape of the 
9(<p) curves. 

3.3. Tilted Halo Potential 

Dekel and Shlosman (1983) and Toomre (1983) 
proposed that observed warps of galaxies may result 
from the fact that the dark matter halo is oblate and 
is rigidly tilted with respect to the inner disk. To 



consider this possibility we generalize equation (55) 
to the case of N tilted rings in a halo potential which 
is rigidly tilted by an angle, say, 6/j = 9^ x (with 
9\ x <C 1) with respect to the z— axis of the inner disk. 
We find 



dt 2 



'dt 



N 



Al J i j (e j -Q h )-Y i c jk (Qj-e k 



where A 2 d - 



Q 2 dj and A 



(61) 



hj 



"hzj 



-Q 2 hj . The 



damping terms are still given by equations (56) and 
(57). 

For a special disk tilt, O^ = funct.(j), the right- 
hand-side of equation (61) is zero so that O^ is time- 
independent. This Qj 1 corresponds to the Laplacian 
surface of the disk (BT, p. 413). It is given by 



N 

£ 



M jk ef; = A 2 hj Q h , 



M 



jk 



= ( A% + K j+ 



it ° jk 



-Cjk(l — 6jk) , 
(62) 

where 5jk is the Kronccker delta and the prime on 
the summation means that the diagonal terms are 
omitted. The damping terms (56) and (57) are of 
course zero for the Laplacian tilt. If the self- gravity 
of the rings (Cjfc) is negligible, 6^ = [A 2 hj /{A 2 hj + 
Adj)]@h- Equation (62) can always be inverted to give 
Qj if A 2 > 0. This is because the N— ring system is 
stable for A 2 > (see discussion following equation 
(53b)) which implies det(Mjk) > 0. 

A general solution of equation (61) can evidently 
be written as Qj(t) — Qf + &j{t), where 6^ obeys 
equation (55) which has no reference to the halo tilt. 
Over a sufficiently long period of time, damping given 
by equation (56) and/or (57) will give Q'j(t) — > as 
discussed in the paragraph after next. 

Figure 19a shows the Laplacian tilt Qj for a repre- 
sentative case with and without the self-gravity of the 
rings. Figure 19b shows the Briggs plots 9(ip) at times 
t = 4 and 8 Gyr for a disk started from a Laplacian 
tilt and from a 20% deviation from a Laplacian tilt. 
The Briggs plot for the Laplacian tilt is in contrast 
with the leading spiral wave observed, for example, in 
M 83 (Figures 1 and 2) (Briggs 1990). However, the 
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Laplacian tilt may be relevant to cases such as NGC 
3718 which shows a relatively straight line-of-nodes. 

There is no evident reason for a disk to be ini- 
tially 'set up' in a Laplacian tilt. Figure 19b shows 
that deviations from the Laplacian tilt evolve to give 
a line-of-nodes which is not straight. As mentioned, 
dissipativc torques due to dynamical friction [equa- 
tion (56)] and/or relative friction [equation(57)] act to 
damp out the deviation (0^) from the Laplacian tilt 
over a period of time. For example, for /3j/Oi = 0.1 in 
equation (56) and no relative friction, the maximum 
line-of-nodes deviation (^49) at 8 Gyr is reduced by a 
factor ~ 0.4 from the case with no dissipation shown 
in Figure 19b. On the other hand for (3'j/Cli = 1.0 and 
no dynamical friction, the maximum line-of-nodes de- 
viation is reduced by a factor ~ 0.5. 

3.4. Initially Tilted Disk Plane 

Here we consider the possibility that at some intial 
time an outer gaseous disk is formed which is tilted 
with respect to the plane of the inner disk. This sit- 
uation could arise by the capture, tidal-disruption, 
radial-spreading, and cooling of a low mass gas cloud 
by a disk galaxy. For a spherical cloud of mass M c i 
and radius a c i in an approximately circular orbit of 
radius r, tidal breakup occurs roughly for a c i >^ 
f {GM cl r 2 /v 2 )i « 3kpc(M c/ /lO lo M )*(r/2Okpc)§ 
(200km/s/w o )3 , where v is the circular velocity in 
the flat rotation region of the galaxy [equation (6b)]. 
Treatment of the cloud disruption is beyond the scope 
of the present work. Instead, we study the initial 
value problem where at t = there is a smooth tran- 
sition from an untilted inner disk to a tilted outer disk 
with a straight line-of-nodes. Specifically, we take 

Qj(t = 0) = ^ jl - cos[ 7 r(r j - n)/(r N - n)]| , 

(63a) 

d®j{t = 0) 



dt 



= 



(636) 



where 9 xo is the initial tilt angle of the outer disk. 
Similar behavior is found for other smooth variations 
of 6j . The smooth variation of 6j avoids the excita- 
tion of short radial wavelength modes. 

Figure 20 shows the nature of the warp evolution 
resulting from the initial conditions (63). Figure 21 
shows a surface plot of the disk at t = 2 Gyr. A 
leading spiral wave appears in the Briggs plot (panel 
c), and its amplitude grows with time. Progressively, 
the interior region of the disk (say, 10 < r < 25 kpc) 



flattens to the plane of the inner disk while the exte- 
rior region warp amplitude 6 grows. The dependences 
shown in Figure 20 for t = 2 — 4 Gyr are qualitatively 
similar to those observed for M 83 (see Figure 2). 
The outer part of the disk shows an inverse power 
law dependence of the line-of-nodes pnr (panel b 
of Figure 20) which is qualitatively similar to that for 
M 83 where cp oc r -1,5 . 

The dissipative torques due to dynamical friction 
[equation (56)] and that due to relative friction [equa- 
tion (57)] affect the results of Figure 20 in different 
ways. For a dynamical friction coefficient /3j/Oi = 
0.1, and no relative friction, the maximum warp am- 
plitude at t — 4 Gyr is reduced by about 30%. On the 
other hand, for no dynamical friction, but a relative 
friction coefficient {3'j = 1.0, the warp at t = 4 Gyr is 
essentially unchanged. 

4. Continuum Limit 

In the continuum limit, the Hunter and Toomre 
(1969) equation for the vertical displacement of the 
disk /i(r, 0, i) can be obtained from equation (55) 
by recalling that h = r9 x (t)sm(j) — r9 y (t)cos(f) — 
Xm[r O exp(i(/>)]. We omit for the moment the damp- 
ing terms (56) and (57). Multiplying equation (55) by 
r exp(i(j)) and taking the imaginary part gives 



h s = - ij (n 2 + a 2 ) h 



N N 

- X! C i k h i + X! C i k r M r k-,^, t)/r k , (64) 

where the left hand side of this equation follows from 
the same steps as for equation (21), where hj = 
h(rj,(f),t), and where Ij = 2irr^ ArY,(rj) with Ar 
taken as a constant for simplicity. The continuum 
limit is TV — > 00 and Ar oc \/N — ► 0. 

From the definitions (14) and (18), we have on the 
right hand side of equation (64), fl 2 + A 2 = R 2 Z + Q, 2 — 
fl 2 . From the decomposition of the potential (§1), it 
is clear that VL 2 — Cl 2 is the contribution to the disk 
rotation from the outer disk; that is, 



n 2 -n 2 = 



ld<P od 



r dr 

E(r') [r - r'cos(^)] 
[r 2 + (r') 2 - 2rr'cos{<j> l )Y 



(65) 
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where d 2 r' = r'dr'dcp', and where the integration is 
over the outer part of the disk, r > r\ , with n chosen 
to be significantly larger than the scale-length of 
the inner disk (see §2). Note that the summations in 
equation (64) are also over the outer part of the disk. 
Thus, one term which contributes to the right hand 
side of equation (64) is 

-I^-^) hj = - Gljhj [ d 2 r> g f--^)] 

J r [. . .J 2 

(66a) 

where E' = S(r'), and the factor in the denominator 
to the 3/2 power is the same as in equation (65). 
Using the definition of Cjk in equation (34), the first 
sum on the right hand side of equation (64) can be 
converted to an integral in the continuum limit, 



E, 



7- 



r \. . . 2 



(666) 



Similarly, the continuum limit of the second sum in 
equation (64) gives 



J2 =GI,jd 2 r' r h{r '-" J] ms( ' "» 

= Glj [ d 2 r' 

J [r 2 + (r') 2 — 2rr'cos 



r 2 + (r') 2 — 2rr'cos((f) — </>')]? 
E' h(r',(j)',t) 



(66c) 

in view of equation (31c). 

Combining equations (64)- (66) and dropping the j 
subscripts gives 



h = -HRi h 



GE / d 



E' (h - ti) 



W , (67) 

[r 2 + (r 1 ) 2 — 2rr'cos(0 — (j)')] 2 

where h' = h(r',<p',t). Note that R z includes the 
vertical restoring force of the inner disk and the 
halo. Thus, equation (67) is the same as Hunter and 
Toomre's (1969) equation for h, in that they did not 
separate out the inner disk and did not include a halo 
potential. 

The Newtonian drag term (56) due to dynamical 
friction gives a contribution 



(68) 



on the right hand side of equation (67). The damping 
term (57) due to relative friction between the rings 



gives a contribution 
1 d 



dr 



dr \r dt J 



(69) 



on the right hand side of equation (67). 

Multiplying equation (67) by 2d 2 r(dh/dt), inte- 
grating over the outer disk, and following the steps 
of HT (their Appendix B) gives, in the absence of 
dissipation, the constant of the motion 



■^warp 



d 2 r E 



^] +A 2 h 2 



G 



drr 2 Y, dr' (r'f E 



J(r,r') 



h(r',<j>,t) h(r,4>,t) 



(70) 



where J(r, r') is given by equation (30b). Substituting 
h(r,(p,t) — r8(t)sm[(f) — (p(t)] from equation (8a), one 
readily finds that £ war p is identical to £ r i ngs = 7i of 
equation (53b). Note that E warp is non-negative and 
therefore the disk is stable if A 2 (r) > 0. 

Including the dissipation terms (68) and (69) gives 



dE, 



warp 



dt 



r v t E 



dr\r dt J 



(71) 



The relative friction is of course zero if h/r is a func- 
tion only of time which corresponds to a rigid tilt. 

The second constant of the motion V rings, found 
in §2.1 and §3 [equation (53c)], can of course be ex- 
pressed in terms of h{r, <f>, t) in the continuum limit. 
Instead, we obtain this new constant directly from 
the HT equation for h by multiplying equation (67) 
by 2d 2 r(dh/d<j>) and integrating over the outer disk. 
This gives 



^Pwarp — 2 / d r E 



. dhy dhdh 



I 



= d 2 rZ r 2 6 2 (n-<p) 



Const. 



(72) 



and V W arp — V rings in the continuum limit, in the 
absence of dissipation. With dissipation included, we 
have not found a simple expression for dV war p/dt. 
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5. Conclusions 

This paper develops a representation for the anti- 
symmetric small amplitude warp dynamics of a self- 
gravitating disk in terms of N tilted concentric rings. 
That is, we consider only azimuthal mode number 
of the warp m = 1 (or —1). This representation is 
suggested by the kinematic ring model of Rogstad 
et al. (1974) which is used in interpreting HI disk 
warps (Briggs 1990). The rings are considered to be 
in a fixed oblate (or prolate) halo potential and in the 
potential of an inner, untiltcd disk. Different initial 
value problems are studied using our N— ring model. 

We first consider in detail the tilting dynamics of 
one ring in the potential of the halo and inner disk. 
The equations of motion are shown to have a par- 
ticularly simple form when written in terms of the 
complex tilt angle 9 = #exp(— iip), where 9 is the ac- 
tual tilt angle, and ip is the angle of the line-of-nodes 
relative to the inertial axis. The single ring has a 
slow and a fast precession mode which are analogous 
to the normal modes of vibration of a non-rotating 
mechancial system. For the slow (fast) mode 0(t) ro- 
tates in the clockwise (counter clockwise) direction in 
the complex plane. The dynamical equation for of 
course has an energy constant of the motion, £ r %ng- 
Additionally, we show that there is a second constant 
of the motion. By writing out the Lagrangian for one 
ring, C = C(6,0,<f), this constant of the motion is 
recognized as the canonical angular momentum for 
the ip coordinate, P v = dC/dip = Const. 

We examine the influence of the drag on the motion 
of one ring through the halo matter due to dynami- 
cal friction. The full calculation of this drag (Nelson 
and Tremaine 1995) is beyond the scope of this work. 
Instead, we include a linear drag torque in the equa- 
tion of motion for which we estimate following the 
approach of Weinberg (1985). Inclusion of this term 
shows that dissipation can destabilize the ring tilt- 
ing in a prolate halo potential. The instability is of 
the negative energy type; that is, it occurs only for 
< 0. 

We next consider the particularly interesting case 
of two tilted, gravitationally interacting rings of dif- 
ferent radii in the potential of the halo and inner disk. 
We derive dynamical equations for Oi = #iexp(— ifi) 
and 62 = 9 2 cxp(— if 2)- This case shows the phe- 
nomenon of phase-locking of the line-of-nodes of the 
two rings. The gravitational interaction between the 
two rings is appropriately measured by the coupling 



strength £12 of equation (48). As £12 increases the 
ring motion changes from that of independent pre- 
cession for £12 < 1 to phase-locked precession for 
£12 > 1 where the lines-of-nodes coincide approxi- 
mately, fx w ipi. This phase-locking persists for the 
case of N > 2 rings. 

We examine the influence of dissipative forces on 
the tilting motion of two rings. In addition to dynam- 
ical friction, which is included as linear drag torques 
on each ring, we include the 'relative' friction which 
can result from the differential vertical motions of two 
gaseous rings due to turbulent viscosity. The corre- 
sponding drag torques on each ring have a definite 
form owing to the requirement that the total angular 
momentum of the two rings be conserved in the ab- 
sence of external torques. The turbulent viscosity is 
assumed to be described by the a model of Shakura 
(1973) and Shakura and Sunyaev (1973). We find that 
the relative friction, as well as the dynamical friction, 
can destabilize the motion of two rings in a prolate 
halo potential. 

We study the case of two tilted rings which rotate 
in opposite directions. This situation, which is per- 
tinent to observed counter-rotating galaxies (see, for 
example, Jore, Broeils, and Haynes 1996), is found 
to be unstable in a prolate halo potential for certain 
conditions in the absence of dissipation. That is, a 
dynamical instability may occur. 

For the general case of ./V— tilted, gravitationally 
interacting rings of radii rj = n + (j — l)Ar in the 
potential of the halo and inner disk we obtain equatins 
of motion for 0j = 9jexp(— ifj) for j = l,..,N. 
The Lagrangian has the form L = £(6j,ifj,6j,tpj), 
and this gives, in the absence of dissipative torques, 
both an energy constant of the motion and a total 
canonical angular momentum constant of the motion, 
P v = ^jdC/ 'difj = Const. We comment on the nu- 
merical solutions for 0j for different values of N, and 
we give a sample power spectrum for N = 49 which 
shows the distribution of the 2N mode frequencies. 
We argue that for treatment of intial value problems 
of time duration t max , mode frequency differences 5oj 
smaller than l/t max are irrelevant. This leads to the 
estimate N = O(10 2 ) for the number of rings needed 
to treat warps in galactic disks. 

In the continuum limit N — > 00, the ring model 
is shown to give the Hunter and Toomre (1969) dy- 
namical equation for the vertical displacement of the 
disk h{r,(j),t). Dissipative torques due to dynamical 
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friction and/or relative friction (for gaseous rings) are 
shown to cause the energy of the warp £ warp to de- 
crease with time. 

We have numerically solved the dynamical equa- 
tions for Qj for N = 49 for four different types of ini- 
tial conditions which may give rise to observed warps 
in galaxies: (1) warp excitation by a passing satellite 
with relatively large impact parameter; (2) excitation 
by a slowly sinking compact minor satellite; (3) warp 
evolution in a tilted halo potential (Dekel and Shlos- 
man 1983; and Toomre 1983); and (4) warp evolution 
resulting from an initially tilted disk plane due to the 
tidal breakup of a gas cloud. The nature of the disk 
response at different times is most clearly shown in 
the polar plots of 6j versus (pj which we refer to as 
Briggs plots. This is one of the forms used by Briggs 
(1990) to describe the warp geometry as deduced by 
fitting a kinematic ring model (Rogstad et al. 1974) 
to HI observations. 

For case (1) we find that the polar plots 0j(<fj) 
have the shape of a leading spiral wave qualitatively 
of the form observed, but the warp amplitude 0j is 
smaller than observed in many galaxies. The small 
warp amplitudes found here are in agreement with 
the conclusions of HT. The leading spiral in the 0j (<pj ) 
curves in the outer part of the disk (where the self- 
gravity between the rings is small) results from the 
dominance of the slow precession mode which gives 
very roughly ipj ~ u> P j cx — 1/r, where u P j is the single 
ring slow precession frequency. The fast precession 
modes disappear in a few Gyrs due to phase-mixing. 

For case (2) of a sinking satellite, we find that the 
0j(ifj) curves have a large noise component and are 
unlike the observed curves in the absence of relative 
friction between the rings. This is due to the violent 
excitation from the satellite passing through the disk. 
However, with relative friction included, the 6j(<fij) 
curves are smooth leading spiral waves qualitatively 
similar in shape to the observed curves. The warp am- 
plitude, although significantly larger than for case (1), 
is still smaller than that for many observed warps. We 
find that the inner part of the disk where 9j is small 
has an approximately straight line-of-nodes (ipj inde- 
pendent of j) for both cases (1) and (2) as observed 
by Briggs (1990) (his 'rule of behavior' No. 2). We 
show that this is due to the above-mentioned phase- 
locking of the lines-of-nodes of the inner rings of the 
disk due to self gravity between the rings. 

For case (3) of an oblate halo potential rigidly tilted 
with respect to the inner disk of the galaxy, there is 



a unique time independent Laplacian warp Qj = Gj" 
(BT, p. 413) which has a straight line-of-nodes (ipj = 
Const.). However, it is unlikely that a galactic HI 
disk is 'set up' with the Laplacian tilt. Deviations 
from Oj evolve to give a line-of-nodes which is not 
straight. Over the age of the galaxy, the deviations 
from may damp out due to dynamical friction 
and/or relative friction between the rings. This case 
is possibly relevant to the warp of NGC 3718 which 
shows an approximately straight line-of-nodes out to 
r w 1.75R Hol (Briggs 1990). 

Case (4) of an initially tilted disk plane with a 
straight line-of-nodes gives: (1) polar plots 9j{<fj) for 
t ~ 2 — 6 Gyr with the leading spiral shape qualita- 
tively of the form observed; (2) straight line-of-nodes 
(ifj w Const.) in the inner part of the disk as ob- 
served as a result of phase-locking due to self-gravity; 
and (3) a warp amplitude linearly dependent on the 
initial tilt of the disk plane. Further, the line-of-nodes 
angle <p(r) in the outer part of the disk is found to 
have an inverse power law dependence on r qualita- 
tively of the form of shown for M 83 (see Figure 2b). 
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Fig. 1. — Surface plot of the warped disk of M 83 
(NGC 5236) obtained from the data of Briggs (1990). 
The rotation of the galaxy is clockwise so that the 
warp of the disk is a leading spiral wave (Briggs 1990). 
The scales are in units of 10 kpc for a Holmberg radius 
of M 83 of RhoI = 7.3 arcmin and a distance of M 83 
of 5.9 Mpc which assumes a Hubble constant H Q = 
65 km/s/Mpc (Briggs 1996, private communication). 
The colors are in 16 steps uniformly spaced between 
Zmin and z max . 

Fig. 2. — Nature of the warp of the disk of M 83 from 
Briggs (1990). The top panel (a) shows the radial (r) 
dependence of the tilt angle 0(r). The middle panel 
(b) shows the r-dependence of the angle of the line-of- 
nodes (p(r). The bottom panel (c) shows the depen- 
dence of 9 on <p, which we refer to as a Briggs plot. 
The radii are in units of 10 kpc for the parameters 
mentioned in the caption of Figure 1 . In these units 
the Holmberg radius is RhoI ~ 1-26. In panel (b), the 
dashed curve shows the least-square-fit of a power law 
dependence to the curve for 1 < r < 3.4 which gives 
ip w 216°r- 1 - 48 (with an R value of 0.993). Note that 
in the Briggs plot ip has been taken to increase in 
the clockwise direction to account for the fact that 
the galaxy rotates in the clockwise direction. Except 
for our dicussion of counter rotating rings, we assume 
that the matter of the rings (representing the galaxy 
disk) rotates in the counter-clockwise direction. Note 
that the data points given here have only an approx- 
imate correspondence with those given by Briggs in 
that we have scanned and digitized the data from his 
figures. 

Fig. 3. — Geometry of one tilted ring of a disk galaxy. 
For the case shown, the line-of-nodes is along the 
x— axis, that is, ip = 0. Also, 9 X > 0, 6 y = 0, and 
n denotes the normal to the ring plane. 

Fig. 4. — Illustrative single ring orbits for a case 
shown in (a) where the slow mode is dominant (C\ = 
0.2, C*2 = 0.8) and in (b) where the fast mode is 
dominant (C\ — 0.8, C 2 = 0.2). For both plots, i 
corresponds to the starting point, and / to the end 
point at t f = 62.8/0, and A 2 /ft 2 = 0.2. 

Fig. 5. — Dependence of the dimensionless integral 
J' 12 defined in equation (45c) on S = r\/r<2. The 
dashed lines correspond to approximations discussed 
in the text. 

Fig. 6. — The top panel (a) shows the dependences of 



the two slow precession frequencies of two co-rotating 
rings in an oblate halo on the strength of the grav- 
itational interaction measured by uo g i2 obtained the 
full equation (43) and the approximate equation (44). 
The frequencies are measured in units of fti. For 
this figure, fti/ft 2 = 2, h/h = 1, and A 2 = 0.2ft 2 . 
[Note that the two fast precession modes, 0^3 and 
L04, have frequencies 1.05 and 2.1 for u> g i2 = 0, and 
1.19 and 2.19 for uo g i2 = 0.1. The uj 3 mode has 
62/61 3> 1 and thus involves mainly motion of ring 
2, whereas the UJ4 mode has —62/61 <C 1 and in- 
volves motion mainly of ring 1.] 

— The bottom panel (b) shows the dependences of 
the amplitude ratios 62/61 for the two slow preces- 
sion modes. For small values of 0J g i2, the u)\ mode in- 
volves motions mainly of ring 2, whereas the u>2 mode 
involves motion mainly of ring 1. 

Fig. 7. — Illustrative orbits {©i(i), 6 2 (i)} for two co- 
rotating rings in an oblate halo for weak (a, top panel) 
and strong (b, bottom panel) coupling for the same 
parameters as Figure 6. The initial conditions are 
6i/6 2 = 1/2 and 6i = = 6 2 . The points labeled 
by f's all correspond to the final time tf = 307r/fti 
(which is arbitrary). The constants of the motion 
£rings and P r ings are constant in the numerical inte- 
grations with fractional errors < 10~ 6 . For the top 
panel, w s i 2 /fti = 0.005 (£12 ~ 0.3), and the rings 
precess almost independently. For the bottom panel, 
w 9 i 2 /fti ~ 0.015 (£12 ~ 0.9), and the rings precess to- 
gether with their lines-of-nodes phase-locked as shown 
in Figure 8. The dividing value separating indepen- 
dent and phase-locked precession is w g i 2 /fti ~ 0.010 
or £12 ~ 0.60. Gf course, this value depends on the 
inital value 6i/6 2 for 6,(0) = 0. 

Fig. 8. — Azimuthal angle difference of the line-of- 
nodes of two rings ip 2 — ipi for cases (a) and (b) of 
Figure 7. At the final time tf = 307r/fti, for the 
phase-locked case (y> 2 — </?i))/ ~ —24°, whereas for 
the unlocked case (<p 2 — tpi)f w 253.5°. 

Fig. 9. — The top panel (a) shows the dependences of 
the two slow precession frequencies of two co-rotating 
rings in an prolate halo on the strength of the gravita- 
tional interaction measured by oj g i2 obtained equation 
(43) . The frequencies are measured in units of fti . For 
this figure, fii/fi 2 = 2, h/h = 1, and A 2 = -0.2ft 2 . 

— The bottom panel (b) shows the dependences of 
the amplitude ratios 62/61 for the two slow preces- 
sion modes. For small values of w g i 2 , the u)\ mode 
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involves motions mainly of ring 2, whereas the 0J2 
mode involves motion mainly of ring 1. For large u> g \2 
the u>2 mode frequency approaches the limit given in 
equation (47) and 62/61 — > 1. 

Fig. 10. — Region of instability of two rings, one 
co-rotating > 0) and the other counter-rotating 
(f^2 < 0), in a prolate halo potential as discussed be- 
low equation (51). 

Fig. 11. — The top panel (a) shows the dependences 
of the two slow precession frequencies of two counter 
rotating rings (f2i > and Q 2 < 0) in an oblate 
halo on the strength of the gravitational interaction 
measured by co g i2 obtained equation (43). The fre- 
quencies are measured in units of fli . For this figure, 
Oi/Q 2 - -2, /1//2 = 1, and A 2 = 0.2Q?. 

— The bottom panel (b) shows the dependences of 
the amplitude ratios 6 2 /6i for the two slow preces- 
sion modes. 

Fig. 12. — The top panel (a) shows the dependence 
of the three slow precession frequencies of a three ring 
system on the strength of the gravitational interaction 
measured by 0J g i2- The frequencies are measured in 
units of For this figure, (ri,r 2 ,r3) = (1,2,3), 

(fii,n 2 ,n 3 ) = (1,5,3), h/h - 1, h/h = 1, 

and A 2 = 0.2S1 2 . For small values of u) g i2, the 
lu\ mode involves motion mainly of ring 3, the oj 2 
mode motion of ring 2, and the mode motion 
of ring 1. For larger values of uo g \ 2 the motions of 
the different rings become coupled. For large uo g i 2 
the slowest mode (oji) approaches the limiting value 
given by the generalization of equation (47) to three 
rings. [The three fast precession modes have fre- 
quencies 0.698, 1.05, and 2.09 for uo g i 2 = and 
0.842, 1.30, and 2.20 for w g i 2 /fii = 0.1.] 

— The bottom panel (b) shows the r-dependence of 
the vertical displacements h (along a line _L to the 
line-of-nodes) for the three slow precession modes for 
0J g i2/&>i = 0.1. The vertical scale is arbitrary except 
for the condition h 2 <C r 2 . Smooth lines have been 
drawn through the calculated values indicated by cir- 
cles. 

Fig. 13. — Illustrative orbits Qj(t) for three co- 
rotating rings in an oblate halo for weak (a, top panel) 
and strong (b, bottom panel) coupling for the same 
parameters as Figure 12. The initial conditions are 
6j oc Tj and Qj = 0. The points labeled by f's all 
correspond to the final time tf = 457r/f2 1 (which is ar- 



bitrary). For the top panel, u> g i 2 /Cli = 0.005 and the 
rings precess almost independently. For the bottom 
panel, oj gl2 /£!i = 0.03 and the rings precess together 
with their lines-of-nodes phase-locked as shown in fig- 
ure 14. 

Fig. 14. — Azimuthal angles <fj(t) of the line-of- 
nodes for three rings for cases of: (a) weak {u g \2 = 
0.005Oi), (b) intermediate (uj g i 2 = 0.015Oi), and (c) 
strong coupling (uj g i 2 = 0.03^i). 

Fig. 15. — The figure shows the power spectrum as 
a function of angular frequency ui of 6 X 2s{t) - the 
x— component tilt angle of the middle ring at r = 25 
kpc - for a N — 49 ring disk given a random ini- 
tial perturbation. The frequency u is normalized by 
the angular velocity of the disk at n, Q\. The spec- 
trum is obtained using a 1024 point FFT. The initial 
perturbation is taken to have a 'white noise' radial 
wavenumber spectrum; that is, all radial wavenum- 
ber modes are excited. The &j(t) are obtained us- 
ing a code which solves equations (55) with the rings 
equally spaced between n = 10 kpc and — 40 kpc. 
Dissipation is neglected (/3j = = f3'j). The mass of 
the inner disk in equation (4) is Md = 6x 1O 1O M , 
and the disk radial scale is = 4 kpc. The halo 
potential is given by equation (6a) with cllipticity 
e = 0.0871, core radius r a = 15 kpc, and circular 
velocity v a = 200 km/s. The neutral hydrogen has a 
total mass of 1.36 x 10 10 .Mq and is distributed accord- 
ing to equation (5) with Rh = 20 kpc. The separation 
between low and high frequency modes is roughly at 
u)/Cti — 0.75 The lowest frequency of the spectrum 
is w 0.0368 while the next two higher fre- 

quencies are at w 0.0521 and 0.0639. Note that the 
limiting value of lo\ from the generalization of equa- 
tion (47) is lo\ — > Lj^pj/ J2Lj ~ — 0.120fii which 
is significantly larger in magnitude than the u)\ ob- 
served. 

Fig. 16. — The top panel (a) shows a polar plot (a 
Briggs plot) of the tilt angle of the disk 9 in degrees 
as a function of the angle of the line-of-nodes (p for 
the case of a retrograde passage of a compact satellite 
for times t=l, 2, and 3 Gyr after the time of closest 
approach. The plot shows a leading spiral wave qual- 
itatively of the form observed, but with a relatively 
small amplitude compared with many observed warps 
(Briggs 1990). 

— The bottom panel (b) shows a comparison of the 
disk response with and without gravitational interac- 
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tion (Cjfc) between the rings. This figure shows the 
locking of the phase ip in the inner part of the disk. 

— The results shown in (a) and (b) are obtained 
from a code which solves equations (55) including the 
torques of equation (59) with N = 49 rings equally 
spaced between n = 10 kpc and r 49 = 40 kpc. 
Thus the separation between rings is Ar = 0.625 kpc. 
The torques are evaluated numerically at each time 
step. The Newtonian drag is neglected, /3j — in 
equation (56). A small relative friction is included, 
(j'./Q.-t = 0.02 in equation (57). [For a disk half- 
thickness Az = 250 pc and sound speed c s = 10 
km/s, this /3j corresponds to a value of Shakura's 
(1973) viscosity parameter a « 0.0624. This friction 
acts to smooth out sharp 'corners' which exist in the 
f3'j = curve 0((p). Owing to this (3j, Brings decreases 
by 16% between t = 1 and 3 Gyr.] The mass of the 
satellite is 2 x 1O 1O M , its core radius is a = 1 kpc, 
and at closest approach it is at a distance r so = 30 
kpc where it has a speed of 400 km/s and is located 
at an angle 30° above the (x,y) plane. The mass of 
the inner disk in equation (4) is Md = 6 x 10 10 Mq, 
and the disk radial scale is = 4 kpc. The halo 
potential is given by equation (6a) with ellipticity 
e = 0.0871, core radius r = 15 kpc, and circular 
velocity v a = 200 km/s so that (k\ z - n£)/fi£ = 0.2 
for r 2 3> r 2 . The mass of halo matter inside of r = 30 
kpc is 7.4 x 10 10 M Q . The neutral hydrogen has a total 
mass of 1.36 x 10 10 M© and is distributed according 
to equation (5) with Rh = 20 kpc. 

Fig. 17. — The top panel (a) shows a polar plot of 
the tilt angle of the disk in degrees as a function 
of the angle of the line-of-nodes ip for the case of a 
prograde passage of a compact satellite at times t = 1 
and t = 4 Gyr after the time of closest approach. The 
other conditions are the same as given in the caption 
to Figure 16. The plot shows that initially (t = 1 
Gyr) a trailing spiral wave is formed. Later, (t = 4 
Gyr) the trailing wave 'unwraps' and gives rise to a 
leading spiral wave. Analysis of the evolution of the 
tilt excited by the prograde satellite shows that it is 
initially a fast outward propagating wave with radial 
phase velocity v^. h w 24 km/s and radial wavelength 
A r w 11 kpc for r ~ 25 kpc. Roughly, A r oc r 2 5 . 

— The bottom panel (b) shows a comparison of the 
orbits [025x(t), &25y(t)] (in degrees) of ring No. 25 at 
radius r 2 5 = 25 kpc for t = — 1 Gyr for the cases of 
retrograde and prograde satellite encounters for the 
conditions of Figures (15a) and (16a). The circles 
and squares are equally spaced at 0.1 Gyr intervals. 



The (i) indicates t = and the (f's) correspond to 
t = 1 Gyr. For the retrograde case, mainly the slow 
precession mode of the ring is excited, whereas for 
the prograde case mainly the fast precession mode is 
initially excited. 

Fig. 18. — The top panel (a) shows a polar plot of 
the tilt angle of the disk 8 in degrees as a function of 
the angle of the line-of-nodes ip for the case of a ret- 
rograde sinking of a compact minor satellite of mass 
M s = 10 10 M Q at times t = 2 and t = 4 Gyr after 
the sinking (r s = 0). Initially, the rings are unper- 
turbed and the satellite is 'turned on' at a distance 
r s = 50 kpc in the (x,z) plane 30° above the (x,y) 
plane. The in-spiral of the satellite is described by 
equation (60a) with In A = 3. The other conditions 
are the same as described in the caption of Figure 16 
except that the relative friction is larger, f}'j/il\ — 0.2. 
The Newtonian drag is neglected. This value of (3'j has 
the effect of damping out short radial wavelength fea- 
tures of the tilt while not appreciably changing the 
overall 9(ip) curve. 

— The bottom panel (b) shows the Briggs plot for a 
prograde sinking for conditions otherwise the same as 
in (a) above. 

Fig. 19. — The top panel (a) shows the Laplacian tilt 
angle of the rings 9^ = 9j x with and without self- 
gravity of the rings (Cjk) for the case where the halo 
potential [equation (6a)] is rigidly tilted by an angle 
Oh = 9hx = 30° with respect to the z— axis of the 
inner disk of the galaxy. The disk and halo are the 
same as in Figure 16. The are obtained by solving 
equation (62). 

— The bottom panel (b) shows the Briggs plots for 
the Laplacian tilt and that for a small initial de- 
viation from a Laplacian tilt at times t = 4 and 
8 Gyr. The initial deviation is taken to be Qj = 
9f x {l + 0.2sin[ 7 r(r j - n)/2(r N - n)]}. This deviation 
vanishes at rj = r\ and it has a maximum at the outer 
edge of the disk. Dissipative torques are neglected. 
This deviation from 0j" gives rise to a leading spiral 
feature at the outer edge of the disk at t = 8 Gyr. For 
a deviation with the +0.2 factor replaced by —0.2 the 
line-of-nodes angle at the outer edge of the disk ((^49) 
at t = 8 Gyr has a similar magnitude to that in Figure 
19b, but it is negative. 

Fig. 20. — The top panel (a) shows the radial depen- 
dence of the tilt angle 9(r) for different times after 
the initial time {t = 0) when the disk tilt is give by 
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equations (63) with outer disk tilt angle 9 xo = 30°. 
The conditions of the disk and halo are otherwise the 
same as for Figure 16. 

—The middle panel (b) shows the radial dependence 
of the line-of- nodes (p(r). The constant value of tp(r) 
in the inner part of the disk results from phase-locking 
due to the self-gravity of the rings as discussed in §2.3 
and §2.5. The dashed curves correspond to least- 
square fits of a power law to ip(r). For t = 2 Gyr, 
the fit gives <p = 498°r" 1 - 42 for r = 23 - 40 kpc (with 
an R value of 0.96). For t = 4 Gyr, the fit gives 
ip = 1680°r- 1 - 79 for r = 28 - 40 kpc (with an R value 
of 0.97). Note that for M 83 (Figure 2b), the least- 
square fit gives ip = 216°r~ 1 - 48 . 

— The bottom panel (c) shows the Briggs plots for a 
sequence of times. A leading spiral wave forms, and 
its amplitude grows with time. The 6(ip) spiral for 
t = 2 — A Gyr is qualitatively similar to that observed 
for M 83 (see Figure 2c). From results not shown 
here the similarity continues up to t ~ 6 Gyr after 
which time 9 and d9/dr become large at the outer 
edge of the disk. The present theory assumes small 
tilting amplitudes S 2 « 1 for the linearization of the 
equations of motion. The linearity allows an arbitrary 
rescaling Qj — > KQj with K = Const, (or 9j — > K6j 
and ipj — ► (fij) in Figure 20. 

Fig. 21. — Surface plot of the warped disk of Figure 
20 at t — 2 Gyr. To facilitate comparison with Figure 
1 for M 83, the rotation of the disk has been taken 
to be clockwise in this figure. The warp is a leading 
spiral wave the same as for M 83. The colors are in 
16 steps uniformly spaced between z min and z max . 
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